Engineering of arbitrary U(A^) transformations by quantum Householder reflections 
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We propose a simple physical implementation of the quantum Householder reflection (QHR) 
M(ii) — I — 2\v) {v\ in a quantum system of A'^ degenerate states (forming a qunit) coupled simulta- 
neously to an ancillary (excited) state by A'^ resonant or nearly resonant pulsed external fields. We 
also introduce the generalized QHR M(?;;(/p) = I + (^e"^ — l) \v) {v\, which can be produced in the 
same A^'-pod system when the fields are appropriately detuned from resonance with the excited state. 
We use these two operators as building blocks in constructing arbitrary preselected unitary trans- 
formations. We show that the most general V{N) transformation can be factorized (and thereby 
produced) by either A'' — 1 standard QHRs and an A'^-dimensional phase gate, ot N — 1 generalized 
QHRs and a one-dimensional phase gate. Viewed mathematically, these QHR factorizations pro- 
vide parametrizations of the U(A') group. As an example, we propose a recipe for constructing the 
quantum Fourier transform (QFT) by at most A'^ interaction steps. For example, QFT requires a 
single QHR for A^ = 2, and only two QHRs ior N = 3 and 4. 

PACS numbers: 03.67.Lx; 32.80.Bx; 03.67.Dd; 03.67.Hk 



I. INTRODUCTION 

Coherent control of quantum dynamics traditionally 
involves scenarios for transfer of population, complete 
or partial, from one bound initial energy state to an- 
other, single or superposition state, or a continuum of 
states. Such techniques are well developed, particularly 
for two-state and three-state systems, e.g. tt pulses 1|, 
adiabatic passage using one or more level crossings [31, 
or stimulated Raman adiabatic passage (STIRAP) and 
its extensions Q. Essentially all these techniques start 
from a single initial state; such a state can be prepared 
experimentally, e.g. by optical pumping. 

In the same time, in contemporary quantum physics 
implementations of specific propagators are often de- 
manded, for example, some fields in quantum information 
lean heavily on the quantum Fourier transform An- 
other example is quantum state engineering when a sys- 
tem starts in a coherent superposition of states; then one 
must construct the entire propagator, while the above 
techniques provide only some transition probabilities. 

The implementation of such propagators is well under- 
stood and used for qubits, i.e. two-state quantum sys- 
tems, upon which the theory of quantum information is 
primarily built [J]. On the other hand, qunits - A^-state 
quantum systems - offer some advantages. For exam- 
ple, a qubit can encode two continuous parameters: the 
population ratio of the two qubit states and the relative 
phase of their amplitudes. A qunit in a pure state can en- 
code 2(A^ — 1) parameters {N — 1 populations and A^ — 1 
relative phases), i.e. by using qunits information can be 
encoded in significantly fewer particles than with qubits. 
This is beneficial for storing quantum information, which 
can be particularly important if the number of particles 
that can be used is restricted, e.g., due to decoherence 
Furthermore, there are indications that using qunits can 
improve error thresholds in fault tolerant computation. 



Physical realizations of qunit operations in the existing 
proposals [H], however, are difficult to implement. These 
implementations use sequences of U(2) operations, i.e. 
transformations acting at each instance of time upon only 
two of the A^ states of the qunit. The general U(A^) trans- 
formation of a qunit requires 0{N'^) such U(2) operations 
5] ; hence the complexity increases rapidly with the qunit 
dimension A^, which makes qunit manipulations challeng- 
ing, even for qutrits {N = 3). 

In this paper, we show that a general '^{N) transforma- 
tion can be implemented physically in a quantum system 
with only N interaction steps. For this purpose we in- 
troduce a compact quantum implementation, in a single 
interaction step, of the Householder reflection The 
latter is a powerful and numerically very robust unitary 
transformation, which has many applications in classi- 
cal data analysis, e.g., in solving systems of linear alge- 
braic equations, finding eigenvalues of high-dimensional 
matrices, least-square optimization, QR decomposition, 
etc. 0. The Householder transformation, acting upon 
an arbitrary A/-dimensional matrix, produces an upper 
(or lower) triangular matrix by A^ — 1 operations. When 
the initial matrix is unitary, the resulting final matrix 
is diagonal, i.e. a phase gate or a unit matrix. We use 
this propery to decompose an arbitrary U(A^) matrix into 
Householder matrices and hence, design a recipe for phys- 
ical realization of a general U(A^) transformation. 

The quantum Householder reflection (QHR) consists of 
a single interaction step involving N simultaneous pulsed 
fields. In contrast to the existing U(2) realizations of 
qunit transformations, here each Householder reflection 
acts simultaneously upon many states: A'^ states in the 
first step, N — 1 states in the second, etc. This allows 
us to greatly reduce the number of physical steps, from 
0{N'^) in U(2) realizations to only 0{N) in our proposal. 

We introduce two types of QHRs: standard QHR and 
generalized QHR; the latter involves an additional phase 
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factor. The physical reahzations of both use simultane- 
ous pulses of precise areas in a system with an iV-pod 
linkage pattern, the difference being that the standard 
QHR operates on exact resonance, whereas the general- 
ized QHR requires specific detunings. Any unitary ma- 
trix can be decomposed into — 1 standard QHRs and a 
phase gate, or into N generalized QHRs, without a phase 
gate. This advantage of the generalized- QHR implemen- 
tation derives from the additional phase in each step, 
which delivers N additional phases in the end, thereby 
making the phase gate unnecessary. 

This paper is organized as follows. In Sec. |TT] we de- 
fine the standard and generalized QHR gates and propose 
physical implementations. In Sec. [nilwe describe the de- 
compositions of a general U(iV) matrix by means of stan- 
dard and generalized QHRs, which provide the routes for 
realization of an arbitrary V{N) transformation. In Sec. 
?? we apply these decompositions to quantum Fourier 
transforms. The conclusions are summarized in Sec. ??. 




|1) |2) |3> \N) 

FIG. 1: Physical realization of the quantum Householder re- 
flection: TV degenerate (in RWA sense) ground states, forming 
the qunit, coherently coupled via a common excited state by 
pulsed external fields of the same time dependence and the 
same detuning, but possibly different amplitudes and phases. 



II. QUANTUM HOUSEHOLDER REFLECTION 
(QHR) 

A. Standart QHR 

An A^-dimensional quantum Householder reflection 
(QHR) is defined as the operator 



M{v) = I-2\v) 



(1) 



where \v) is an A^-dimensional normalized complex 
column-vector and I is the identity operator. The QHR 
^ is hermitean and unitary, M(w) — M(w)^ = M(z;)"^, 
which means that M(w) is involutary, M^(w) = I; in ad- 
dition, det M(w) = — 1. If the vector \v) is real, M(v) has 
a simple geometric interpretation: reflection with respect 
to an (A'^ — l)-dimensional plane with a normal vector \v) ; 
in the complex case the interpretation is more involved. 
In general, the Householder vector \v) is complex, which 
implies that it contains 2{N — 1) real parameters (taking 
into account the normalization condition and the unim- 
portant global phase). 



B. Generalized QHR 

We define the generalized QHR as 

M{v-^)^l+{e'^ -l)\v){v\, (2) 

where \v) is again an A^-dimensional normalized complex 
column- vector and ip is an arbitrary phase. The standard 
QHR ([T]) is a special case of the generalized QHR ([2]) 
for (p — tt: M(w;7r) = M(i;). The generalized QHR is 
unitary. 



M{v; ipy^ = M{v; ip)' = M{v; -ip), 



(3) 



C. Physical implementations 

1. Coherently driven N-pod system 

The standard and generalized QHRs have simple phys- 
ical realizations. Consider the (A^ -I- l)-state system 
with N degenerate [in the rotating-wave approximation 
(RWA) sense [H] ground states |n) (n = 1,2,...,A^), 
which represent the qunit, coupled coherently and simul- 
taneously by N external fields to an ancillary excited 
state |e) = |A^ -I- 1), as shown in Fig. [1] [8]. Such an A^- 
pod system can be formed, e.g., by coupling the magnetic 
sublevels of several J = 1 levels to a single J = level 
by polarized laser pulses for a qutrit only one J = 1 
level suffices. The propagator UN+i{t,tQ) of this system 
obeys the Schrodinger equation. 



th-VN+lit,to) = ll{t)lJN+l{t,to), (4) 

at 



with the RWA Hamiltonian [H 



ni{t) 

^2 {t) 





n\{t) nut) 



n^it) 
n%{t) 2A{t) 



, (5) 



and the initial condition U^v+i {to, to) = I. The excited 
state |e) can be generally off resonance by a detuning 
A(t) Q, which, however, must be the same for all fields. 
The functions r2i(<), . . . , r2jv(i) are the Rabi frequencies 
of the couplings between the ground states and the ex- 
cited state; we require that they have the same time de- 
pendence, described by the envelope function / (t) , but 
we allow for different amplitudes Xn ^^'^ phases 



and its determinant is det M 



njt)^Xnf{t)e'^- {n = l,2,...,N). 



(6) 
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By using the Morris-Shore transformation |9| the cou- 
pled {N + l)-state system can be decomposed into a 
set of TV — 1 dark ground states, which are superposi- 
tions of qunit states, and a two-state system, consist- 
ing of a bright ground state and the excited state |e) 
0. This two-state system is driven by a Hamiltonian 



involving the same detuning A{t) as in Eq. (O, and the 
coupling is the root-mean-square (rms) Rabi frequency 

The exact solution for the propagator reads [1] 



U 



N+l — 



XiX2e 



3l2 



XiXjye 



(a-1) 

_5*Xi£ 



(a-1) 



X2Xn±_ 



hi 



(q _ 1) xix«f ^" 5^ 



,X2e 



Hi 



(7) 



l + (a-l)^ h 



Here x = \J'Ylin=iX^ is the rms peak Rabi frequency 
and — Pk~ 0m {k,m — 1,2,..., TV) are the relative 
phases of the external fields. The complex parameters 
a and h (with = 1 — |a|^) are the Cayley-Klein pa- 
rameters of the SU(2) propagator for the Morris- Shore 
bright-excited two-state system. 



2. Standart QHR: exact resonance 

In the case of exact resonance (A — 0) the Cayley- 
Klein parameters for any pulse shape f{t) are 



A 



'% sm ■ 



(8) 



iV-dimensional QHR vector \v) are the normalized Rabi 
frequencies, with the accompanying phases. 



1 

X 



(12) 



Hence the propagator U'^ within the degenerate iV-state 
manifold of the iV-pod system driven by the Hamiltonian 
^ , with A = and rms pulse area © , represents indeed 
a physical realization of QHR in a single interaction step. 
Any QHR vector p2|) can be produced by appropriately 
selecting the peak couplings x„ and the phases while 
obeyng Eq. pH)) (e.g., by adjusting the pulse duration). 



3. Generalized QHR 



where A is the rms pulse area. 



If 



A^X f{t)dt. 



A = 2 (2fc 1) TT (fc = 0,1,2, 



(9) 



(10) 



then a — —1, sin (A/2) = 0, and the last row and col- 
umn of the propagator ([7]) vanish, except for the diagonal 
element, which is —1; the propagator ([7]) reduces to 



U 



JV+l 



r ^0 
U'^ : 

L J 

••• 0-1 



(11) 



Here U'^ is an A^-dimensional unitary matrix (with 
det U'^ = —1), which represents the propagator within 
the A^-state degenerate manifold; it has exactly the QHR 
form (P), U'' = M.{v\ tt) = M{v). The components of the 



The unitary propagator ([7]) for a = e*"^ (|6| = 0) re- 
duces to 







••• e-''^ 



(13) 



where, as is easily verified, we have U"^ = M(w; ip), and 
hence, the propagator U"^ represents a physical realiza- 
tion of the generalized QHR ([2]). The vector \v) is again 
given by Eq. (fT2|) . The condition a = e^'^ for ip ^ 0,t: 
can only be realized off resonance (A ^ 0). There is a 
beautiful off-resonance solution to the Schrodinger equa- 
tion - the Rosen-Zener (RZ) model - which we shall use 
here to exemplify the generalized QHR. 

The Rozen-Zener (RZ) model [13] can be seen as an ex- 
tension of the resonance solution ^ to nonzero detuning 
for a special pulse shape (hyperbolic-secant), 



fit) 
A(t) 



sech(t/T) , 
Ao. 



(14a) 
(14b) 
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The Cayley-Klein parameter a reads [1, [T3| 



r (i + + i*Aor) r (i - + i*AoT) ' 

(15) 

where r(z) is Euler's gamma function. Using the reflec- 
tion formula F (i + z) F (i — z) = tt/ costtz, we find 



sin^ (iTTxr) 
cosh^ (ittAqT) ■ 



(16) 



Hence in this model, \a\ = 1 for xT = 21 
{I = 0, 1, 2, . . .); then the last row and the last column of 
the propagator ([7]) vanish, except the diagonal element. 



The phase ip of a 



depends on the detuning Aq and 



for an arbitrary integer I we find from Eq. (115 



and hence 



n 

fc=0 



Aor + z(2fc+l) 
AoT-i(2fc+l) 



2 arg [| [AoT + z (2fc + 1)] 



(17) 



(18) 



fe=0 



This can be seen as an algebraic equation for Aq, which 
has / real solutions. For example, for ^ = 1 [which cor- 
responds to rms pulse area A — 27r], we have AqT — 
cot((/3/2). Hence the generalized-QHR phase (p can be 
produced by an appropriate choice of the detuning Aq. 

The use of nonresonant interaction, besides providing 
an additional phase parameter, has another important 
advantage over resonant pulses: lower transient popu- 
lation of the intermediate state. This can be crucial if 
the lifetime of this state is short compared to the interac- 
tion duration. Equation (I18|l provides the opportunity to 
control this transient population, which is proportional 
to A-2, by using large peak Rabi frequency (implying 
larger /) and find the largest solution for A. It is im- 
portant that the standard QHR can also be realized off 
resonance, by selecting a detuning Aq for which Lp = tt. 



III. QHR DECOMPOSITION OF V{N) 

A. Standard-QHR decomposition 

We shall show that QHR is a very efficient tool for 
constructing a general V{N) qunit gate. In particular, 
we shall show that any A^-dimensional unitary matrix U 
(U~^ = U^) can be expressed as a product of iV — 1 
standard QHRs M(t)„) (n = 1, 2, — 1) and a phase 
gate * {<j}i,<j}2, ■ ■ • ,0Ar), 

U = M{VI)M{V2) ■ ■ ■ M{VN-1)^ 0w) , (19) 

where 

* ((/-i , 02, . . . , 0^) = diag(e^^i , e^<^= , . . . , e'^« ) . (20) 



We shall prove this assertion by explicitly construct- 
ing the decomposition The standard QHRs M(i;„) 
involve vectors |wn), which we construct as follows. First 
we define the normalized vector as 



"1 



^2 [l-Re(uiie-»^i)]' 



(21) 



where the vector \u„) denotes the nth column of U = 
{ukn}, 4>i = arguii, and |ei) = [1,0, ...,0] . We find 

M{vi)\u^) = e'^Mei), (22a) 
m{vi)\u.^) = \u^)+2e-"'^^uin\vi), (22b) 
(ei|M(wi)|u„) =0 (n = 2,3,...,7V). (22c) 

Hence the action of M(wi) upon U nullifies the first row 
and the first column except for the first element. 



M(ui)U 



e«0i 
^ 











Ujv-1 

L J 



(23) 



where Utv-i is a \]{N — 1) matrix. We repeat the same 
procedure on ]V[(wi)U and construct the vector |i;2), 



J<t'2 



\V2) = 



162) 



^2 [l-Re(u^2e-»^2)]' 



(24) 



where the vector [wa) is the second column of M(?;i)U, 
02 = arg[M(t;i)U]22, and [ea) = [0, 1, 0, . . . , 0]^. The 
corresponding QHR M(i;2), apphed to M(wi)U, has the 
following effects: (i) nullifies the second row and the sec- 
ond column of M(wi)U except for the diagonal element, 
which becomes e*'^^ ^ and (ii) does not change the first row 
and the first column. By repeating the same procedure 
— 1 times, we construct A'' — 1 consecutive Householder 
reflections, which nullify all off-diagonal elements, to pro- 
duce a diagonal matrix comprising N phase factors, 



M(i;Ar_i)...M(«i)U- *((/)!, 



(25) 



which completes the proof of Eq. (fT9|) since M(i;) = 
M.{v)-^. If U is a SU(iV) matrix then det * = ±1, 

meaning X^^Li 0„ = or tt. 

We note that the choice of the QHRs M(z;„) is not 
unique; for example, the first QHR M(i;i) can be con- 
structed from the first row of U, instead of the first col- 
umn. Furthermore, the final diagonal matrix (I20p occurs 
due to the unitarity of U, which leads to Eq. (|22cp : a 
QHR sequence produces a triangular matrix in general. 

The QHR decomposition (fT9|) of the U(iV) group into 
iV — 1 Householder matrices ([T]) and a phase gate provides 
a simple and efficient physical realization of a general 
transformation of a qunit by only TV — 1 interaction steps 
and a phase gate; this is a significant advance compared 
to 0{N^) operations in existing recipes. Each QHR vec- 
tor is iV-dimensional, but the nonzero elements decrease 
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from N in to just 2 in |'i;7v_i), and so does the num- 
ber of fields required for each QHR, see Eq. p^ . 

The decomposition (fT9|) is also of mathematical inter- 
est because it provides a very natural parametrization of 
the V{N) group. Indeed, a QHR vector with n nonzero 
elements contains 2{n — 1) real parameters (because of 
the normalization and the irrelevant global phase). The 
phase gate ([20)) contains N phases. Hence Eq. (ITO)) in- 
volves J2n=2 2(n— 1)+A^ = N"^ real parameters, as should 
be the case for a general V{N) matrix. 



B. Generalized-QHR decomposition 

We shall show now that any unitary matrix U can 
be expressed as a product of N generalized QHRs 
M {vn, </?n) (n = 1, 2, . . . , A^) defined by Eq. (O, without 
a phase gate, that is 



N 



u = n M{v„ 



(26) 



We first define the normalized vector 
, , 1 /2sin((^i/2) 



^^1 



1 



l-ui 



«i)-|ei)), (27) 



where the vector |u„) denotes again the nth column of 
U and ipi = 2 arg (1 — un) — tt. It is readily seen that 

Mivr,-^i)\ui) = |ei), (28a) 
{ei\Mivi;~ip^)\ur,) = (n = 2, 3, . . . , TV ).(28b) 

Therefore, the action of M(wi; —tpi) upon U nullifies the 
first row and the first column except for the first element, 
which is turned into unity. 



M(vi;-(pi)U 



1 

^ 





(29) 



where Ujv_i is a U(iV — 1) matrix. We repeat the same 
procedure on Uat-i and construct the vector 



1 



^^2 



g-JVa _ 1 



^ 2 Sin ((^2/2) 

|l-"22l 



(K)-|e2)), (30) 



where the vector 



is the second column of 



M(t;i; —ipi)XJ and — 2 arg (1 — 1*22) — tt. The action 
of M(w2; — V2) upon M(wi; —(pi)XJ has the following ef- 
fects: (i) nullifies the second row and the second column 
of 'M.{vi; —(pi)XJ except for the diagonal element which 
is turned into unity, and (ii) does not change the first 
row and the first column of M(wi; —ipi)\J. By repeating 



the same procedure N times, we construct N consecu- 
tive generalized Householder reflections, which nullify all 
off-diagonal elements to produce the identity matrix. 



n M(f„ 



)U 



(31) 



By recalling Eq. ([3]) we obtain Eq. immediately. 
Note that the last QHR M.{vn] ^m) — *(0, • ■ • , 0, ip^) is 
actually a one-dimensional phase gate. 

Therefore the use of generalized QHRs replaces the 
iV-dimensional phase gate needed in the standard-QHR 
implementation l|19p by a one-dimensional phase gate 
$(0, . . . , 0, (/3^). We point out that again, as for the 
standard QHRs M(w„), the choice of any of the gen- 
eralized QHRs M(w„;(/5„) is not unique because it can 
be constructed from the respective row, rather than the 
column, of the corresponding matrix. 



C. Examples 



1. Quhit 

As an example of the QHR decomposition we first 
consider the qubit, which is the conventional system for 
quantum information processing. The conventional re- 
alization of a general U(2) transformation involves three 
interactions: two phase gates and one rotation R (i?) 



(32) 



U = *(ai,a2)R(i?)*(0,a3) 



Already for a qubit, the QHR implementations ([T 
and ((26)) arc superior to Eq. ((32)) because they only re- 
quire one QHR and one phase gate. 



U = M(«)*(0i,(/)2), 
U = M(^;;^i)* (0,(^2) 



(33a) 
(33b) 



2. Qutrit 

As a second example we consider a qutrit — a three- 
state quantum system. The most general transforma- 
tion of a qutrit belongs to the U(3) group, which can be 
parametrized by nine real parameters; respectively, the 
SU(3) group is described by eight real parameters. A 
SU(2) factorization of SU(3) reads [HI 

U = R23 (ai,/3i,7i)Ri2 (a2,/32,a2)R23 {^3,^3,73) > 

(34) 

where are SU(2) subgroups of SU(3), with the SU(2) 
submatix occupying the mth and nth rows and columns 
of Rmn- Hence this implementation ((M)) of SU(3) re- 
quires three SU(2) gates, each involving three qubit gates 
((32)). i.e. nine qubit gates in total (which can be reduced 
to seven by combining adjacent phase gates). With the 
present QHR implementation ((55)) of SU(2) the number 
of operations can be reduced to six. 
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Already for SU(3) or U(3), the present QHR imple- 
mentations and (pS)) are considerably more efficient 
because they require only two QHRs and a phase gate, 

U = M(t;i)M(t;2)*(^i,(/.2,03), (35a) 
U - M(z;i;(^i)M(z;2;¥'2)*(0, 0,(^3). (35b) 

As an example, the arbitrarily chosen SU(3) gate 



U = 



0.864e~2'^*/3 0.282ei5"/i9 0.4166"^''*/^ 
0.382e°i'"'^* 0.902e^'^*/" 0.203e°-*°8^'' 
0.327e-o'^897r» o.328e4'^'/5 0.886e°"35^'' 



(36) 

(keeping 3 significant digits) can be realized with two 
standard QHRs and a phase gate, with 



\vi) = [0.260e'"/^0.734e"■l^°"^0.628e-•'■^«9"']|^37a) 
\V2) = [0,0.651e-°-i34-,0.759e*^-"i°"]^, 

,-0.6677ri „0.8667ri „-0.1997ril 



* = diag{e-''-667-,e°-«66x\e-°-i39«|_ 



(37b) 
(37c) 



Alternatively, the same SU(3) gate can be realized 
by two generalized QHRs and a phase gate (j35bp . with 
ifii = — 0.6937r, = 0.6537r, 1^93 = 0.047r, and 



[0.955e°-3°^",0.226e 
[0,0.987e''-347'^\0.161e- 



-0.7077ri 



,0.193e' 



0.3647ri 



-0.3837ril 



(38b) 



B. Qutrit 

For a qutrit the QFT matrix reads 



1 1 1 

-j^ p27ri/3 g — 27ri/3 
-j^ g — 27r'i/3 g27r'i/3 



1 - \/3,l,l 



l^'2; = 



'1 + V2 
2V2 



0, 1 - V2, 



The ^enerafeerf-QHR decomposition reads 
Uf = M(«i;7r)M(z;2;7r/2), 



1 



vi 



V2 



•1 



1 
71 



1- V3, 1,1 



1 r iT 



(42) 



The standard-QHR decomposition reads 

Uf = M(vi)M(v2)*(0,7r/4,-37r/4) (43a) 



(43b) 
(43c) 

(44a) 
(44b) 

(44c) 



Here the first QHR M(ui;7r) = M(ui) is the same for 
the standard- and generalized- QHR implementations. 



IV. QUANTUM FOURIER TRANSFORM 

The quantum Fourier transform (QFT) is a key in- 
gredient in quantum factoring, quantum search, general- 
ized phase estimation, the hidden subgroup problem, and 
many other quantum algorithms [^] . The QFT is defined 
as the unitary operator with the following action on an 
orthonormal set of states \n) (n = 1, 2 . . . , N): 



1 ^ 

fe=i 



A. Qubit 

For a qubit, TJ-'^ is the Hadamard gate 



1 



1 1 

1 -1 



(40) 



which can be written as a single QHR, = M(u), with 

(41) 



1 



-V2-V2,V2 + V2 . 
Here the standard and generalized QHRs coincide. 



C. Quartit 



For a quartit {N = 4) the QFT matrix reads 



U5 



1111 

1 i -1 -i 
1-11-1 
1 -i -1 i 



(45) 



The standard-QHR decomposition reads 

M(z;i)]V[(?;2)*(0, 7r/4, 0, -37r/4), (46a) 

(46b) 

T 



l^^l) = i [-1,1,1, if, 



\V2) = 



2V2 



0,1- \/2,0, 



The generalized-QHR decomposition reads 

Uf = M(,;i;7r)M(«2;7r/2), 
\vi) = ^[-1,1,1, if, 

\V2) = -^[0,1,0,-lf . 



(46c) 

(47a) 
(47b) 

(47c) 



Again, the ffist QHR M(wi;7r) = M(wi) is the same 
for the standard- and generalized- QHR implementations. 
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FIG. 2: Deviation T,f,k=i |(f^Jv),fe - {U^)j^\ of tlie prop- 
agator U]v(f) from the QFT matrix versus time for 
A'' = 2, 3, 4, for generalized- QHR implementations. The 
pulses for N = 2 are centered at time r = —ST, whereas 
ioT N = 3 and 4 at times r = — 5T and r = ST. We have 
assumed sech pulse shapes (|14a|l and rms pulse area A = 2n 
(x = 2). The individual couplings Xn given by the com- 
ponents of the generahzed-QHR vectors (|4T|) for N = 2, (|44b|) 
and (fiic)) for = 3, ()47b|l and (|i7c)) for iV = 4, each multi- 
phed by x- AH phases /3,j are zero. The detunings are A = 
for the first steps and A — 1/T for the second. 



Interestingly, the QFT for = 4 is decomposed with 
only two QHRs, rather than three, without phase gates. 

Figure ?? shows the time evolution of the propaga- 
tor U7v(0 towards the respective QFT matrix U^, for 
N — 2,3,4, for realizations with generalized QHRs. As 
time progresses, the deviation of XJN{t) from van- 
ishes steadily in all cases. As predicted, QFT is realized 
with just a single QHR for N = 2 and with just two 
QHRs for iV = 3 and 4. 



V. DISCUSSION AND CONCLUSIONS 

We have proposed a simple physical implementation 
of the quantum Householder reflection in a coherently 
driven A^-pod system. We have shown that the most gen- 
eral U(A'^) transformation of a qunit can be constructed 
by at most iV — 1 standard QHRs and an A^-dimensional 
phase gate, or by A^ — 1 generalized QHRs (each hav- 
ing an extra phase parameter compared to the standard 
QHR) and a one-dimensional phase gate, i.e. by only N 
physical operations. This significant improvement over 
the existing setups [involving 0{N'^) operations] can be 
crucial in making quantum state engineering and opera- 
tions with qunits experimentally feasible. 

The Householder gate is superior already for a qubit 
because the general U(2) gate needs just two gates, a 
QHR and a phase gate, compared to three gates in exist- 
ing implementations. For a qutrit, the QHR realization of 
U(3) requires only three gates, compared to at least seven 
hitherto. The QHR implementation of the U(iV) gate is 



particularly important for qutrits because of the straight- 
forward physical implementation inaJ=l<->J = 
transition (Fig. [T|); the results, of course, apply to any 
A'^, and can be accomplished, for instance, by using more 
J = 1 levels. 

We have given examples for QHR implementations of 
quantum Fourier transforms. The QHR realization of 
QFT for a qubit requires a single interaction step, com- 
pared to two steps hitherto. The QHR realization is par- 
ticularly efficient for a quartit (A^ — 4), where the QFT 
is synthesized with only two QHR gates [as for a qutrit 
[N = 3)], much fewer than 0(4^) in the existing SU(2) 
proposals. The components of the Householder vectors 
are the amplitudes of the respective couplings. It is im- 
portant that all QHR phases are relative phases of the 
external control fields, e.g. relative laser phases, which 
are much easier to control than dynamic and geometric 
phases. 

The generalized QHR requires off-resonant pulsed in- 
teractions, appropriately detuned from resonance. The 
standard QHR can be realized both on and off resonance. 
The off-resonance implementation has the advantage that 
only negligible transient population is placed into the 
(possibly decaying) ancillary excited state; however, it 
requires a specific value of the detuning. 

In the existing SU(2) proposals, each interaction step 
involves a single SU(2) (or Givens) rotation. The differ- 
ence between the Givens rotation and the Householder 
reflection is that, when applied to an arbitrary matrix, 
the Givens rotation nullifies a single matrix element; the 
Householder reflection nullifies an entire row (or column). 
When the matrix is unitary, a single Householder reflec- 
tion nuUifles one column and one row simultaneously. 
Hence the Householder reflection is N times faster than 
the Givens SU(2) transformation. 

In atoms and ions qubits are encoded usually in degen- 
erate ground sublevels, and the coupling between them 
is accomplished by off-resonant interactions, via an inter- 
mediate state, which is eliminated adiabatically to pro- 
duce an effective Raman coupling. In doing so, the phase 
relation between the two Raman fields is lost. In our pro- 
posal we use resonant, or nearly-resonant, fields; no adi- 
abatic elimination is performed and the phase relation is 
preserved in the resulting QHR propagator. Therefore, 
already for N — 2, the QHR contains an additional phase 
parameter compared to previous realizations, which re- 
duces the number of steps for U(2) operations from 3 to 2. 
Hence, even for a qubit there is a clear improvement. It is 
also significant that resonant interactions, which we use, 
require less interaction energy than off-resonant interac- 
tions; this may be crucial in the case of weak couplings. 

We conclude by emphasizing that the wide-spread use 
of the Householder reflection in classical data analysis 
promises that the proposed quantum implementation has 
the potential to become a powerful tool for quantum state 
engineering and quantum information processing. 
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